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Abstract. This paper studies the problem of reconstructing binary ma- 
trices that are only accessible through few evaluations of their discrete 
X-rays. Such question is prominently motivated by the demand in ma- 
terial science for developing a tool for the reconstruction of crystalline 
structures from their images obtained by high-resolution transmission 
electron microscopy. Various approaches have been suggested for solving 
the general problem of reconstructing binary matrices that are given by 
their discrete X-rays in a number of directions, but more work have to 
be done to handle the ill-posedness of the problem. We can tackle this 
ill-posedness by limiting the set of possible solutions, by using appropri- 
ate a priori information, to only those which are reasonably typical of 
the class of matrices which contains the unknown matrix that we wish 
to reconstruct. Mathematically, this information is modelled in terms of 
a class of binary matrices to which the solution must belong. Several 
papers study the problem on classes of binary matrices on which some 
connectivity and convexity constraints are imposed. 

We study the reconstruction problem on some new classes consisting of 
binary matrices with periodicity properties, and we propose a polynomial- 
time algorithm for reconstructing these binary matrices from their or- 
thogonal discrete X-rays. 

keyvirords: combinatorial problem, discrete tomography, binary matrix, 
polyomino, periodic constraint, discrete X-rays. 



1 Introduction 

The present paper studies the possibility of determining the geometrical as- 
pects of a discrete physical structure whose interior is accessible only through a 
small number of measurements of the atoms lying along a fixed set of directions. 
This is the central theme of Discrete Tomography and the principal motivation 
of this study is in the attempt to reconstruct three-dimensional crystals from 
two-dimensional images taken by a transmission electron microscope. The quan- 
titative analysis of these images can be used to determine the number of atoms 



in atomic lines in certain directions |22 2^. The question is to deduce the local 



atomic structure of the crystal from the atomic hne count data. The goal is to 
use the reconstruction technique for quality control in VLSI (Very Large Scale 
Integration) technology. Before showing the results of this paper, we give a brief 
survey of the relevant contributions in Discrete Tomography. 

Clearly, the best known and most important part of the general area of to- 
mography is Computerized Tomography, an invaluable tool in medical diagnosis 
and many other areas including biology, chemistry and material science. Com- 
puterized Tomography is the process of obtaining the density distribution within 
a physical structure from multiple X-rays. More formally, we attempts to recon- 
struct a density function f{x) for x in or from knowledge of its line integral 
Xf{L) = Jj^ f{x)dx for lines L through the space. This Hne integral is the X-ray 
of f{x) along L. The mapping f ^ Xf is known as the Radon transform. The 
mathematics of Computerized Tomography is quite well understood. Appropri- 
ate quadratures [ p8|j29[ | of the Radon inversion formula are used, with concepts 
from calculus and continuous mathematics playing the main role. 

Usually, the physical structure has a very big variety of density values, and 
so a large number of X-rays are necessary to ensure the accurate reconstruction 
of their distribution. In some cases the structure that we want to reconstruct 
has only a small number of possible values. For example, a large number of 
objects encountered in industrial computerized tomography (for the purpose of 
non-destructive testing or reverse engineering) |9| are made of a single homoge- 
nous material. In many of these applications there are strong technical reasons 
why only a few X-rays of the structure can be physically determined. Discrete 
Tomography is the area of Computerized Tomography in which these special 
cases are studied. The name Discrete Tomography is due to Larry Shepp, who 
organized the first meeting devoted to the topics in 1994. 

An example of such a case is the above-mentioned problem of determining local 
atomic structure of a crystal from the atomic line count data. In a simple but 
highly relevant model suggested by Peter Schwander and Larry Shepp the pos- 
sible atom locations in a unit cell of a crystal are defined on the integer lattice 
Z'^, while the electron beams are modeled as lines parallel to a given direction. 
The presence of an atom at a specific location corresponds to a pixel value one 
at the location; the absence of an atom corresponds to a pixel value zero. The 
number of atoms along certain lines through the crystal (i.e., the sum of pixel 
values along those lines) define the discrete X-rays of the atomic structure with 
respect to the family of lines involved. Since in practice, one degree of freedom 
for moving the imaging device is used to control the position of the crystal, the 
view directions for which data are provided lie all in the same plane. This means 
that the 3D-problem leads itself to a 2D-slice-by-slice reconstruction. Therefore, 
the crystal is represented by a binary matrix and its discrete X-ray along a di- 
rection u is an integral vector giving the sum of its elements on each line parallel 
to u (see Fig. ||). 

Measurements are usually only available along two, three or four directions, 
which is much less than what is typical used in Computerized Tomography (a 
few hundred). In fact, the electron microscope makes measurements at the atomic 
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Fig. 1. A subset _F of with the corresponding binary matrix. X^i^o^F, X^q^^F and 
are the discrete X-rays in the directions (1, 0), (0, 1) and (1, 1). 



level and uses high energy (and so deeply penetrating) rays which can corrupt 
the crystal itself. So, it can take only two, three or four images of the crystal 
before the energy of the radiations destroys it, or at least changes permanently 
its atomical configuration so that the subsequent radiations will see something 
different from the original one. 

Now, the problem is to invert the discrete Radon transform, i.e., to reconstruct 
the binary matrix from this small set of discrete X-rays. More precisely, the 
basic question is to determine, given a set of directions ui, . . . ,Uk and a set of 
integral vectors Xi , . . . Xk , whether there exists a binary matrix F whose discrete 
X-rays along ui, . . .uu are Xi, . . . , Xk. The general methods of Computerized 
Tomography cannot be used effectively if the number of X-rays is so small, and 
they seems unlikely to work in practice. 

Discrete Tomography has its own mathematical theory mostly based on dis- 
crete mathematics. It has some strong connection with combinatorics and ge- 
ometry. We wish to point out that the mathematical techniques developed in 
Discrete Tomography have applications in other fields such as: image process- 
ing Q , statistical data security |^ , biplane angiography , graph theory 
and so on. As a survey of the state of the art of Discrete Tomography we can 
suggest the book | |2C| ]. 

Interestingly, mathematicians have been concerned with abstract formulations of 
these problems before the emergence of the practical applications. Many prob- 
lems of Discrete Tomography were first discussed as combinatorial problems 
during the late 1950s and early 1960s. In 1957 Ryser |Q and Gale |l3) gave a 
necessary and sufficient condition for a pair of vectors being the discrete X-rays 
of a binary matrix along horizontal and vertical directions. The discrete X-rays 
in horizontal and vertical directions are equal to row and column sums of the 
matrix. They gave an exact combinatorial characterization of the row and col- 
umn sums that correspond to a binary matrix, and they derived a fast 0{nm) 
time algorithm for reconstructing a matrix, where n and m denote its sizes. We 



refer the reader to an excellent survey on the binary matrices with given row 
and column sums by Brualdi 

The space of solutions of the reconstruction problem, however, is really huge and 
in general quite impossible to control. A good idea may seem to start increasing 
the number of X-rays one by one in order to decrease the number of solutions. 
Unfortunately, the reconstruction problem becomes intractable when the number 
of X-rays is greater than two, as proved in [Q. This means that (unless P = NP) 
exact reconstructions require, in general, an exponential amount of time. In 
polynomial time only approximate solutions can be expected. In this context, 
an approximate solution is close to the optimal one if its discrete X-rays in 
the set of prescribed directions are close to those of the original set. Various 
approaches have been suggested for solving the problem . Recently, an 

interesting method | [T^ for finding an approximate solutions has been proposed. 
Even though the reconstruction problem is intractable, some simple algorithms 
proposed in JlSf have good worst-case bounds and they perform even better in 
computational practice. 

Unluckly, this is not still enough. During the last meeting devoted to Dis- 
crete Tomography, Gabor T. Herman [|l9j and Peter Gritzmann stress the 
fact that various approaches have been suggested for solving the general prob- 
lem of reconstructing binary matrices that are given by their discrete X-rays 
in a small number of directions, but more work has to be done to handle the 
ill-posedness of the problem. In fact, the relevant measure for the quality of a 
solution of the problem would be its deviation from the original matrix. Hence in 
order to establish this deviation we would have to know the real binary matrix. 
However, the goal is to find this unknown original binary matrix so we can only 
consider measures for the quality of a solution based on the given input dis- 
crete X-rays. We have a good solution in this sense if its discrete X-rays in the 
prescribed directions are close to those of the original matrix. Unfortunately, if 
the input data do not uniquely determine the matrix even a solution having the 
given discrete X-rays may be very different from the unknown original matrix. 
It is shown in |^ that extremely small changes in the data may lead to en- 
tirely different solutions. Consequently, the problem is ill-posed, and in a strict 
mathematical setting we are not able to solve this problem and get the correct 
solution. 

In most practical application we have some a priori information about the im- 
ages that have to be reconstructed. So, we can tackle the algorithmic challenges 
induced by the ill-posedness by limiting the class of possible solutions using ap- 
propriate prior information. The reconstruction algorithms can take advantage 
of this further information to reconstruct the binary images. 

A first approach is given in ||2^, where it is posed the hypothesis that the 
binary matrix is a typical member of a class of binary matrices having a certain 
Gibbs distribution. Then, by using this information we can limit the class of 
possible solutions to only those which are close to the given unknown binary 
matrix. A modified Metropolis algorithm based on the known Gibbs prior pro- 



vides a good tool to move the reconstruction process toward the correct solution 
when the discrete X-rays by themselves are not sufficient to find such solution. 

A second approach modelled a priori information in terms of a subclass of 
binary images to which the solution must belong. Several papers study the prob- 
lem on classes of binary matrices having convexity or connectivity properties. 
By using these geometric properties we reduce the class of possible solutions. 
For instance, there is a uniqueness result jisf for the subclass of convex binary 
matrices (i.e., finite subsets F of Z" such that F = Z" n conv{F)). It is proved 
that a convex binary matrix is uniquely determined by its discrete X-rays in 
certain prescribed sets of four directions or in any seven non-parallel copla- 
nar directions. Moreover, there are efficient algorithms for reconstructing binary 
matrices of these subclasses defined by convexity or connectivity properties. For 
example, there are polynomial time algorithms to reconstruct hv-convex polyomi- 
noes (i.e., two-dimensional binary matrices which are 4-connected and 

convex in the horizontal and vertical directions) and convex binary matrices 
from their discrete X-rays. At the moment, several researchers are studying the 
following stability question: given a binary matrix having some connectivity and 
convexity properties and its discrete X-rays along three or four directions, is it 
possible that small changes in the data lead to "dramatic" change in the binary 
matrix ? 

In this paper, we take the second approach into consideration, and we propose 
some new subclasses consisting of binary matrices with periodicity properties. 
The periodicity is a natural constraint and it has not yet been studied in Dis- 
crete Tomography. We provide a polynomial-time algorithm for reconstructing 
(p, l)-periodical binary matrices from their discrete X-rays in the horizontal and 
vertical directions (i.e., row and column sums). The basic idea of the algorithm 
is to determine a polynomial transformation of our reconstruction problem to 2- 
Satisfiability problem which can be solved in linear time |^ . A similar approach 
has been described in P,p^. 

We wish to point out that this paper is only an initial approach to the 
problem of reconstructing binary matrices having periodicity properties from a 
small number of discrete X-rays. There are many open problems on these classes 
of binary matrices of interest to researchers in Discrete Tomography and related 
fields: the problem of uniqueness, the problem of reconstruction from three or 
more X-rays, the problem of reconstructing binary matrices having convexity 
and periodicity properties, and so on. 

2 Definitions and preliminaries 

Notations. Let A™^" be a binary matrix, = '^i-j ^^-d Cj — '^i.j^ ^'^^ 

each 1 < i < m and 1 < j <n. We define R = (ri, . . . , r™) and C = (ci, . . . , c„) 
as the vectors of row and column sums of A, respectively. The enumeration of 
the rows and columns of A starts with row 1 and column 1 which intersect in 
the upper left position of A. A realization of {R, C) is a matrix B whose row 
and column sums are R and C. 



A binary matrix A™^" is said to be [p, q) -periodical if j = 1 implies that 

ai+p_j^q = lifl<i+|5<m and 1 < j + q < n, 

ai-pj-q = 1 if l<i— p<m and 1 < j — q < n. 
Such a matrix is said to have period (p, q). 

For any given couple (x, y) such that a^^y = 1 we define the set P of propagation 
of the value in position {x,y) in direction {p,q), P = {(x + kp,y + kq)\\ < 
X + kp < m,l < y + kq < n,k e Z}. Such set is called a line. Each line has 
a starting point, which is its leftmost point, and an ending point, which is its 
rightmost point. We say that a line starts on column j and ends on column j' 
when its starting and ending points are on column j and j' , respectively. 
The notion of box is a crucial part for our work. Let A be a {p, (7)-periodical 
matrix. From the periodicity it follows that if there exists an index 1 < i < m 
such that 

fi — Ti+j} + k, then the positions on row i, from column n — q-\rl to column n, 

contain at least k elements equal to 1. Such positions form a box at the end 

of row i and will be addressed to as right box (rt); 
Ti + k = ri+p, then on row i+p, from column 1 to column q we have k elements 

equal to 1. Such positions form a box at the beginning of the row i and will 

be addressed to as left box (It); 

We define the upper and lower boxes (up and Iw respectively) on columns in the 
same way (see Fig. ^), if there exists an index I < j < n such that 

Cj — Cj^q + k then the positions on column j, from row m — p + 1 to row m, 
contain at least k elements equal to 1. Such positions form a box at the end 
of column j and will be addressed to as lower box (Iw); 

Cj -\- k = Cj+q then the positions on column j + q, from row 1 to row p, contain 
at least k elements equal to 1. Such positions form a box at the beginning 
of column j and will be addressed to as upper box (up); 
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Fig. 2. A (2, l)-periodical binary matrix with upper (up), lower (Iw), right (rt) and 
left (It) boxes. 



Definitions of polyominoes. A polyomino P is a finite union of elementary cells 
of the lattice Z x Z whose interior is connected. This means that, for any pair 
of cells of P there exists a lattice path in P connecting them (see Fig. |^(a)). A 
lattice path is a path made up of horizontal and vertical unitary steps. These 
sets are well-known combinatorial objects and are called digital 4-connected 
sets in discrete geometry and computer vision. We point out that a polyomino 
can be easily represented by a binary matrix. 

A polyomino is said to be v-convex [h-convex], when its intersection with any 
vertical [horizontal] line is convex. A polyomino is hv- convex or simply convex 
when it is both horizontal and vertical convex. A parallelogram polyomino is a 
polyomino whose boundary consists of two non intersecting paths (except at their 
origin and extremity) having only north or west steps. Fig. ^ shows polyominoes 
having the above-mentioned geometric properties. 



a) 



b) 



c) 



d) 



Fig. 3. a) A polyomino. b) A h,-convex polyomino. c) A hv-convex polyomino. d) A 
parallelogram polyomino. 



3 Periodicity (1,1) 

Let A be a (l,l)-periodical matrix. 

By definition of boxes for p — 1 and q = 1 the boxes are reduced to only one cell 
and the integer k of the definition takes only the values or 1. If there exists an 
index l<i<moTl<j<n such that 

ri = r.i+1 + 1 then Oi^n 1- 
ri + 1 = ri+i then 0^+14 = 1. 
Cj = Cj+i + 1 then a„ij — 1. 
Cj + 1 = Cj+i then aij+i = 1. 

A preprocessing part uses the previous box properties to extract the fixed part 
(called F) of the reconstruction matrix. The following algorithm performed on 
a given pair of vectors (i?, C) gives, if a solution exists, the fixed part (namely 
the matrix F) and a pair of vectors (i?', C") of the mobile part. 

Propagation(x, y, F, £) 



P = {{x + k,y + k)\l < X + k < m,l < y + k < n,k E Z}; 

For all e P do P,,, = l,rf+'' = rf^ - l,cf+'^ = cf^ - 1; 

The main program finds the fixed I's by considering the differences between 
the values of the pair of vectors. For each fixed 1, the procedure Propagation 
fills by periodicity (1,1) the matrix F and decreases the row and column sum of 
the current matrix. 

At the end of the preprocessing part both row and column vectors are homoge- 
neous (r^ = p for all 1 < i < m and = 7 for all I < j < n). 
Now, either (i?', C) are zero vectors and then the solution is unique and equal 
to F, or we perform a reconstruction from homogeneous X-rays with periodicity 
(1,1)- 

Since the vectors are homogeneous we can extend the periodicity on a torus. 
Indeed, suppose that = r.^^^ with = X;"=i a'ij and r'^_^_^ = YTj=i - By 
periodicity a^^ = fli+i j+i, for j — 1, . . . ,n ~ \, implies that aj„ = a^+i^i, for 
i = l,...,m— 1. In other terms the values of the matrix A' are mapped on a 
cylinder. The same argument in column proves that the values of the matrix A' 
are mapped on a torus. That is if a ■ ^ = 1 then a^^^ m.j+i mod n = 1 and 

^2 — 1 mod m,j — l mod n ^' 

So, a solution is formed by loops, namely a beginning in (i,j) with a'^ j — 1 and 
a propagation by periodicity (1,1) until the position All the loops have 

the same length. As the vectors are homogeneous, we can compute the number 
of loops. Using this strong condition and the algorithm of Ryser in the first 
row in order to place the loops, we can reconstruct easily a solution in 0{mn) 
time. 

Another remark is the arithmetical nature of the stability of the solution. We 
can prove that if m and n are relatively prime then there is only one solution. 
Indeed in this case, to perform a reconstruction of a binary matrix with homoge- 
nous vector, we have only one loop felling the whole matrix and then cither the 
matrix A' is full of I's or full of O's and nothing between because of the toric 
conditions. 

Proposition 1. Let R G and C e N". // gcd(n, m) = 1, then there is at 
most a {1,1) -periodical matrix having row and column sums equal to {R,C). 

For example, if we perform a reconstruction with a matrix with m rows and 
m + 1 columns and periodicity (1, 1) then the solution is unique if it exists. 

Example 1. For R — (2, 2, 1, 2), C — (2, 1, 2, 2), the algorithm gives the matrix 
10 
00 1 
00 
1000 

and the vectors R = (1, 1, 1, 1), C = (1, 1, 1, 1). 

We can reconstruct two solutions for R — (2, 2, 1, 2), C — (2, 1, 2, 2): 



10 10 
0101 
00 10 
10 1 

11 
10 1 
10 0' 
10 10 

For R = (1,2,2,2),C = (2,1,2,1,1), the algorithm gives a fixed part and 
R = (0, 0, 0, 0), C = (0, 0, 0, 0, 0), then the solution is unique. 
00 100 
10010 

1 r 

10 100 

4 Periodicity (p, 1) with 1 < p < m 

Let A be a matrix with periodicity (p, 1). 

The preprocessing part uses only the row sums in order to find the fixed part of 
the reconstruction. In fact, by definition of boxes for g = 1 the horizontal boxes 
arc reduced to only one cell and the integer k of the definition takes only the 
values or 1 . If there exists an index 1 <i < m such that 

n = n+p + 1 then Ui^n = 1- 
Ti + l = Ti+i then ai+p,i = 1. 

A preprocessing part uses the previous box properties to extract the fixed 
part (called F) of the reconstruction matrix. The following algorithm performed 
on a given pair of vectors {R, C) gives, if a solution exists, the fixed part (namely 
the matrix F) and a pair of (i?', C") of the mobile part. 

Algorithm 2 

Input: A pair of integral vectors (i?, C); 

Output:If PB = 0, then it gives Matrix F and couple of integral vector (i?', C); 

or If PB = 1 Failure in the reconstruction; 
Initialisation: e := 0, PB := 0, := i?, C^^) ■.= C,F = 0™^"; 
while i?^'^) is positive non homogeneous vector and PB = do 
:= := CW; 

determine first index s.t. r^^^ ^ ''i+pi 

if r^^^-* = r^^^p + 1 then x := i,y := n, Propagation(a;, y, F, £); 

else if rf^ + 1 = r^^^ then x := i + p,y := 1, Propagation (a; , y, F, £) 
else PB := 1; 



and 

A' — 

^2 — 



£:=l+l; 
end while; 



Propagation(a;, y, F, £) 

P = {{x + kp,y + k)\l <x + kp<m,l<y + k<n,keZ}; 
For all e P do Fij = l,rf+'^ = - l,cf+^^ = cf - 1; 

The main program finds the fixed I's by considering the differences between 

the values of the row sums. For each fixed 1. the procedure Propagation fills by 
periodicity {p, 1) the matrix F and decreases the row and column sum of the 
current matrix. 

At the end of the preprocessing part the row vector sum R' of A' has the same 
value on indices in arithmetical progression of rank p: r[ = r^_^_p = r^_|_2p • • ■ = 
where i = L ot L + 1. This set of element of the row sums R' of A' is 
called line of R'. The minimum length of each line of R' is L = [^J . The number 
of lines of length L + 1 and L of R' is n^+i = m mod p and ul = p — ul+i, 
respectively. 

Example 2. If (p, 1) = (2,1), and R = (2, 3, 2, 4, 3, 4, 2), C = (3,4,3,3,4,2,1), 
the algorithm gives the matrix 

0000 100 

0000000 

0000010 
F= 1000000 

1000001 

0100000 

10 

and the new vectors are R' = (1, 3, 1, 3, 1, 3, 1), C = (1, 2, 3, 3, 3, 1, 0). Since m = 
7 and p = 2, we have that L = 3,ni = 1, ul+i = 1. So, R' = (1, 3, 1, 3, 1, 3, 1) 
contains a line of length L = 3 and a line of length L + 1 = 4. Since the lines of 
length L = 3 and L + 1 are (*, 3, *, 3, *, 3, *) and (1, *, 1, *, 1, *, 1), respectively, 
we have that matrix A' contains a line of length three lines of length L = 3 
starting from the second row and L + 1 = 4 starting from the first row. 

Now, we prove now the values of A' are mapped on a cylinder. We have 
= with r- = Y,]=iK,j ^^'^ A+p = YTj=iO-'i+p,j- From the periodicity 
«ij = «i+p,j+i' foi' i = Ir-- - 1, it follows that = <„, for i = 

1,- ■ ■m—p.In other terms the values of the matrix A' are mapped on a cylinder. 

Thus a 1 on the first p rows (at position (x,y),l < x < p,l < y < n) can 
be extended by periodicity on the matrix A' by I's in positions {x + kp,y + k 
mod n) with k = 0, - ■ ■ where £= Lor L + 1. The matrix A' is in particular 
composed on a cylinder of lines in direction (p, 1) of length L or L + 1. In addition 
to that the number of lines of length L + 1 is exactly riL+i = m mod p and the 
number of lines of length L is til = p — til+i- 



4.1 A reduction to the problem of reconstructing a special class of 
h-convex binary matrices lying on a cylinder 

Let A' be a solution for a given {R',C') with r'^ = r'^_^_p — r-+2p ■ ■ ■ — 
where ^ = L or L + 1. We now perform a reduction of reconstruction of {p, 1)- 
periodical matrix A' on a cyhnder to a reconstruction of a special ft,-convex 
matrix A" on a cylinder. 

By the previous construction, matrix A' is formed by lines {x + kp, j + k 
mod n) with k = 0, . . . ,£ — 1 where £ — L oi L + 1. The starting points of the 
lines is the set of position S ~ {{x, y)\l < x < p, 1 < y < n, a'^ y — 1} . 
S is ordered by: (x, y) < {x' , y') if and only if x > x' and y < y' (i.e., we proceed 
from bottom to up and from left to right). 

Let S' be the set 5* with an extra index of the rank in the previous order. Each 
element of S' is a triple {x, y, o), where {x, y) is an element of S and o the rank 
in the order. Now, we can describe the reduction. 

Reduction. Let (x, y, o) G S' . The point (x, y) is the starting point of a line 
of matrix A' having length i or L + 1. This line gives a horizontal bar of I's 
begins in position (o, y + r mod n) with r = 1, . . . , £ where £ = L \i x > n^+i 
or£ = L + lifx< fiL^i. The set of these horizontal bars gives the /i-convex 
matrix A" . 



Notice that, this transformation makes the column sum of A" equal to the col- 



umn sum of A' . In Section 4.3, we will show the inverse reduction that provides 
A' from A" . 

Example 3. Let us take the following matrix A' with periodicity (2, 1) and R' = 
(3, 2, 3, 2, 3) and C = (3, 3, 2, 1, 2, 2) into consideration. 

10 1001 

10 10 
A' = 1 1 1 

1000 1 

110 10 

Since m = 5 and p ^ 2, we have that L = 2, n^+i = 1. The matrix A' is 
composed of three lines of length i + 1 = 3 and two lines of length L = 2. The 
starting points in the first two rows (the two first indices are the position in the 
matrix A' and the last index is the rank in the order) are: 

5' = {(2,1,1),(1,1,2),(1,3,3),(2,5,4),(1,6,5)}. 

The transformation gives the following /i-convex matrix A" mapped on a cylinder 
with three bars of length three and two bars of length two. 

1 1 00 
111000 

A" = 1 1 1 
11 
110 1 

The column sums of A" are equal to C" = (3, 3, 2, 1, 2, 2). 



We point out that the order on the starting points adds the following constraints: 

Condition 1. On each column of A" can start at most bars of length L and 
at most TiL+i bars of length L + 1. Moreover, by proceeding from up to down 
on the column at first we find the bars of the length L and then the bars of the 
length L + 1 (see Fig. ||). 
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Fig. 4. A column of a binary matrix satisfying condition 1. 



We denote the class of /i-convex binary matrices lying on a cylinder and satisfying 
condition 1 by 7iC(nL, n^+i). 

By this property, the matrix A" consists of four disjoint zones B, C, E and P 
whose boundaries are three paths having only north or west steps, and a"(z, j) = 
1 for (i, j) S C U P and a"(j, j) = for (j, j) G BUE (see Fig. | and the matrix 
A" of the previous example). Notice that, the matrices of the class 'HC{nL, Jii+i) 
are set of parallelogram polyominoes lying on a cylinder. 



From the reduction it follows that the problem of reconstructing a {p, 1)- 
periodical binary matrix A' having row and column sums {R',C') (output of 
Algorithm 1) is equivalent to the problem of reconstructing a binary matrix A" 
of 'HV{nL,nL+i) having column sums C", m rows of length L and L + 1. We 
denote this reconstruction problem on the cylinder by RHC problem. 

In the following subsection, we determine a polynomial transformation of 
RHC problem to 2-Satisfiability problem (2-SAT). 



4.2 A reduction to the 2-SAT problem 

Given an instance I of RHC problem, we want to build a 2-SAT formula Q (a 
formula in conjunctive normal form, where each clause has at most two literals) 
whose satisfiability is linked to the existence of a solution for / in such a way: if 
i? is satisfiable, then we are able to reconstruct a solution for I in P-time and, 
vice versa, each solution of / gives an evaluation of the variables satisfying Q in 
P-time. We will do not show the proofs of the lemmas of this section for brevity's 
sake. Let I be an instance of RHC problem; that is: 

- two integers riL, riL+i; 

- a couple (Z/, C), where L and L-\-l are the only possible values of the row sums 
of a binary matrix A" of HC{nL,nL+i), solution of I, and C = (ci, . . . , c„) 
is its column sums. 

- an integer m which denotes the number of rows of A". 

The formula i7 that wc want to construct is the conjunction of three 2-SAT 
formulas: J7i which encodes the geometrical constraints of A", ^2 which gives 
the consistency of A" with the couple {L, C) and, finally, which imposes the 
constraints of condition 1 on each column of A" . The variables of the formula fi 
belong to the union of the four disjoint sots of variables: 

B = {h{i,j) : I <i <m,l < j <n} , C = {c(i, j) : 1 < i < m, 1 < j < n} , 

V = {p{i,j) '■ i < i < in,l < j < n} , £ = {e{i,j) '■l<i<m,l<j< n}. 

We use the variables of the set Xi, with Xi e {B,C,V,£}, to represent the four 
disjoint zones B, C, P, E inside A. 

Coding in fl\ the geometrical constraints of A" . Formula i?i is the con- 
junction of the following sets of clauses: 

r {x{i,j)^x{i-l,j)) A {x{i,j)=^x{i,j + l)) for a; ecu f 
Corners = l\i ^ < 

{ {x{i,j) ^ x{i + A {x{i,j) ^ x{i,j - 1)) foT X G B\JT 

Disj = { iKhi) A ip{ij) => biij)) A (e(i, j) c{i,j)) } 

Compl = /\. .{b{i,j) ^ c{i,j)} 

Anch = {e{l,L) Ap{m,L+l) Ap{r,l)} 
with 1 <r < m. 

Definition 1. Let V\ he an evaluation of the variables in B, C, V, £ which 
satisfies J7i . We define the binary matrix A" of size m x n as follows: 



{c{i,j) = 1 A e{t,j) = 0) ^ a"{i,j) = 1 , p{i,j) = 1 => a"{i,j) = 1 , 
{b{i,j) = 1 A p{i,j) = 0) a"{i,j) = , e{i,j) = ^ a"{i,j) = 0. 



It is immediate to check that A" is weh defined. 

The matrix A" contains the four zones B,C,E and P of yl such that: £ X, 
with X £ {B,C, E, P} if and only if x{i,j) — 1 (see Fig. ||). From Corners, 
Disj, Compl and Anch we deduce the foUowing properties these four zones; 

Lemma 1. i) {B,C} is a partition of A" , P C B and E C C ; 
ii) the boundary of zones B,C,E and P is made up of three paths having only 
north or west steps; 

Hi) there does not exist a column of A" containing both points of P and points 
ofE. 
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Fig. 5. The matrix A" and the four zones B,C,E and P of ^ defined by the formulas 
Corners, Disj, Compl and Anch. 



Coding in f22 the bound of the row and column sums of A". The 

formula is the conjunction of the following sets of clauses: 



LBC - A, 



LBR = A, 



if j > L, e{i,j) b{i + Cj,j) 
if j < L, b{i,j) ^p{i + m-c.j,j) 
if j > L,e{i,j) ^ b{i + Cj,j) 
if j < L, b{i,j) ^p{i + m- Cj,j) 
ifi<r, b{i,i) ^ e{i, j + L + l) 
if i > r, p{i, j) ^ c{i, j + n — L — 1) 
Hi <r, b{i,j) ^e{i,j + L) 
if i > r, p{i,j) c{i,j + n - L) 



The formulas LBC and UBC give a fower and an upper bound for the column 
sums of A. The formula LBR express that the row sums are greater than L. 
Finally, the formula UBR express that the row sums are smaller than L + 1. 
More precisely, 

Lemma 2. Let A" be the binary matrix defined by means of the valuation V2 
which satisfies Qi A o.s in Definition ^. We have that: 

i) the column sums of A" are equal to C ~ (ci, . . . , c„); 
a) each row sum of A" has value L or L + 1 (see Figure 
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L <= L c(i,j) <= L + 1 
j = i 



L <= E p(i,j) + c(i,j) <= L + 1 
j = i 



E p(ij) + c(i,j) = cj E c(i,j) = cj 



Fig. 6. The matrix A" with the bounds on row and column sums. 



Coding in 173 the maximum number of bars of length L and i + 1 
starting on each column of A". The formula is the conjunction of the 
following sets of clauses: 



bbl = A„ <^ 



BBl+1 = A,., < 



if i<m-c„, e(i, j) 6(i - RL+i, j - i - 1) 

if m — c„ < « < r, r — (m — c„) — 1 > n^+i ^ b{m — c„ + 1, n — L) 

if(i>rAj<L), p{i, j) =^h{i - ul+i, j + n - L — I) 

a i < n — Cn, c{i, j) ^ e{i — n^, j + L) 

if n — c„ < z < r, r — {m — Cn) — 1 > zi ^ c{m — Cn + ?^L + l,n — L) 

if i > r, c{i,j)^p{i~nL,j + L-n) 



Lemma 3. Let A" he the binary matrix defined by means of the valuation V3 
which satisfies i7i A i72 A ils as in Definition ^. We have that: 

i) on each column of A" can start at most bars of length L, 

ii) on each column of A" can start at most n^+i bars of length L + 1. 

By Lemmas ^ and ||, matrix A" G HV{nL,nL+i)- By Lemma ||, matrix A" 
satisfies the tomographic constraints. Therefore: 

Theorem 1. f2i A f22 ^ ^3 is satis fiahle if and only if there is a binary matrix 
A" of T-LV{nL^ ^^L+i) having column sum C' , m rows of length L and L + 1. 

Since f2i A X?2 A is a boolean formula in conjunctive normal form with 
at most two literals in each clause, by Theorem Q we have a polynomial time 
transformation of HRC problem to 2-SAT problem which can be solved in linear 
time 



4.3 Final Step 

By performing the previous reduction and an algorithm for solving 2-SAT prob- 
lem we obtain a matrix A" of TCV{nL, til+i) having column sums C", m rows 
of length L and L + 1, where ul+i = rn mod p and = p — n^+i. Now, for 
determining a (p, l)-periodical matrix A' having row and column sums equal to 



(i?', C"), we have to perform the inverse of the reduction defined in Section 4.1 



We point out that (i?',C") is the output of Algorithm 2. This inverse reduc- 
tion should provides A' from A" . The following algorithm describes this inverse 
reduction. 



Algorithm 3 

Input: the matrix A" whose column sums is vector C", and vector R' 

which is homogeneous with respect to p (R' is the output of Algorithm 2); 
Output: the matrix A' having row and column sums {R', C); 

Step 1: determine the two vectors C^^) = {c[^\ ■■■ ,cl^^), C^^) = (cj^\ . . . , d^^) 
such that c^^^ and c|^^ are the number of bars starting from the i-th 
column of A" having length L and L + 1, respectively, ; 

Step 2: construct the matrix A' of size m x n in such a way: 
if 1 < i < (m — p) and 1 < j < we set a'{i,j) = 0; 

if {m— p+l)<i<{m— p + ul+i) and I < j <n, then perform Ryser's 
reconstruction algorithm on the row and column vectors: 

if {m — p + n^+i -I- 1) < J < m and I < j < n, then perform Ryser's 
reconstruction algorithm on the row and column vectors: 

ir,n~p+nL+i+i, 'Tni) and C(2); 
Step 3: for allm — p+1 < i < m and I < j < n, 

if a'{i, j) = 1 then a'{i — pk, (j + k) mod n) = 1 with l<i — pk<m-~p. 



Proposition 2. The matrix A' which is the output of the Algorithm 3 is (p, 1)- 
periodical and it has row and column sums equal to {R',C'). 



We do not show the proof of this Proposition for brevity's sake. 



By performing Algorithm 1, the reduction of Section 4.1, the reduction to 2-SAT 
problem, a linear-time algorithm for solving 2-SAT problem and Algorithm 
3, we obtain a (p, l)-periodical binary matrix A having row and column sums 
equal to {R,C). Since each step can be perfomed in polynomial time, we have 
that: 

Theorem 2. The problem of reconstructing (p, l)-periodical binary matrices from 
their row and column sums can be solved in polynomial time. 



5 Conclusions 

Our main purpose has been to introduce periodicity properties in terms relevant 
for Discrete Tomography. The periodicity is a natural constraint and it has not 
yet been studied in this field. The motivation of this study is in the attempt 
to tackle the ill-posedness of the reconstruction problem by limiting the class 
of possible solutions using appropriate prior information. This means that, we 
modelled a priori information in terms of a subclass of binary images to which 
the solution must belong. 

By using the periodicity properties we reduce the class of possible solutions. 
For instance, we proved a uniqueness result for the class of binary matrices having 
period (1, 1). We have shown a simple greedy algorithm for reconstructing this 
class of matrices from their row and column sums. This reconstruction problem 
becomes more difficult for the binary matrices having period (p, 1) or (1, <?). We 
have described a polynomial-time algorithm for solving this problem which use 
a reduction to 2-Satisfiability problem. We stress the fact that an interesting 
property of this approach is that it can be used for reconstructing parallegram 
polyominoes lying a cylinder from row and column sums. 

The future challenges concern the reconstruction of binary matrices with a 
generical period {p, q) . We wish to point out that this paper is only an initial 
approach to the problem of reconstructing binary matrices having periodicity 
properties from a small number of discrete X-rays. Lot of work should be done 
to understand such environment: we only challenge the reconstruction problem 
from two X-rays in some special cases, but MANY consistency, reconstruction 
and uniqueness problems can be reformulated imposing periodical constraints. 
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